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A generalized model of the spread of the Hantavirus in mice populations is presented on the 
basis of recent observational findings concerning the movement characteristics of the mice that 

■ carry the infection. The factual information behind the generalization is based on mark-recapture 
' observations reported in Giuggioli et al. [Bull. Math. Biol. 67 1135 (2005)] that have necessitated 
I the introduction of home ranges in the simple model of Hantavirus spread presented by Abramson 

^ , and Kenkre [Phys. Rev. E 66 11912 (2002)]. The essential feature of the model presented here is 

O ' the existence of adult mice that remain largely confined to locations near their home ranges, and 

itinerant juvenile mice that are not so confined, and, during their search for their own homes, move 
and infect both other juveniles and adults that they meet during their movement. The model is 
presented at three levels of description: mean field, kinetic and configuration. Results of calculations 
are shown explicitly from the mean field equations and the simulation rules, and are found to agree 
in some respects and to differ in others. The origin of the differences is shown to lie in spatial 

■ correlations. It is indicated how mark-recapture observations in the field may be employed to verify 
' the applicability of the theory. 

O ! PACS numbers: 87.19.Xx, 87.23.Cc, 05.45.-a 

O^! I. BACKGROUND AND MOTIVATION FOR THE STUDY 

' The spread of epidemics is an important topic that has received a great deal of attention from researchers in recent 
^ \ times 0; S 0i IE ■ This interest stems from multiple factors. On the utilitarian side, concerns regarding human 

■ health provide an obvious reason for carrying out such studies. An equally important motivation arises purely from 
intellectual sources: the desire to gain a general understanding of spatially resolved strongly interacting systems on 

I ■ a macroscopic scale. Examples of epidemics of particular interest to the interdisciplinary scientist are the plague, the 
I ' West Nile virus, and the Hantavirus. Prom the many-body aspect in (theoretical) modeling activities, the first presents 
in formidable problems related to the number of carriers and diverse interactions of the fleas and their connections to 
their hosts. The second provides an interesting example of a system vifith two carriers, mosquitoes and birds, with 
Q \ disparate lifespans (weeks and years respectively) and suggests time scale disparity analyses known in condensed 
matter sciences (see e.g. a discussion of the Born-Oppenheimer approximation in solids in ref. |8r). The simplest 
of the epidemics to study, from the conceptual viewpoint of a theorist, is the Hantavirus Q, S S S IE • It is the 
Q^. subject of the study presented below. 

' The details of the Hantavirus epidemic may be found in the review by Yates et al. ^ and related refs. |^,|E^IEI^- 
• w^ , The essential feature of the epidemic from the modeling point of view is that the infection is carried by mice that 
' move from one physical location to another, and is transmitted to other mice through what are probably aggressive 
H ' encounters. It is believed that the mice and the virus have coexisted for millions of years and therefore the mice do 
- - ■ not die, nor are otherwise impaired, from contraction of the virus. This feature is peculiar to the Hantavirus and not 
shared with other epidemics such as the plague or the West Nile virus where the carriers may die from contracting the 
virus. Oscillations in the population which are characteristic of such behavior are therefore absent in the Hantavirus. 
Furthermore, in the Hantavirus context, no mice are born infected: infection may only be contracted from other mice 
after birth: there is no "vertical transmission" of the disease. The human population is incidental to the evolution of 
infection within the mouse population since humans get the virus from the mice but have no feedback effects on the 
mice in the infection process. 

These features of the epidemic led a few years ago to the construction by two of the present authors of a simple 
model of the spread of the Hantavirus in the mouse population. In addition to the above features particular to the 
epidemic, interaction of the mice with the environment through standard logistic terms |^ ITol ] and motion of the mice 
over the terrain viewed as a simple and free random walk were the elements that went into the making of the model. 
The model was introduced at the usual three levels of the description of a many-body system: the most crude, the 
mean field level, the more detailed, the kinetic level, and the most detailed, the configuration level. The last was 
handled through simulations j^. The first two levels were treated in ref. Q, and most active studies were carried 
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out at the kinetic level 's'l which describes through the use of partial differential equations, the time evolution of 
the average but spatially resolved mice populations Mg (susceptible) and Mi (infected). The specific equations, with 
t denoting time and x representing location in a space of appropriate dimensions, a 2-d description being typically 
sufficient, are 

^ . KM. + MO - CM. - "" f; - aM.M. + DV^M., 

ot K(x,t) 

dM, MA Ms + MA , 

— = -cAf,- ; +aM,M, + DV^M,. 1 

ot K{x,t) 

Considerations such as those arising from gender and age of the animals were tacitly neglected, and all parameters 
except K were considered to be independent of time. 

The parameters of the model are a, b, c, K and D. All processes are assumed to be occurring continuously in time, 
an approximation which seems reasonable unless observations or phenomena are particularly chosen to be probed at 
very short time scales. The processes of birth and death of the mice are represented as occurring at rates b and c, 
respectively. The transmission of infection occurs through encounters, the aggression parameter being a. Saturation 
of the population is assumed to occur by competition for the resources among the mice through the environmental 
parameter K which describes nutrition available to the mice. This parameter, whose product with & — c is called 
the 'carrying capacity', controls the so-called logistic term which resolves in the standard manner the paradox of 
Malthusian explosion of the population. 

Let us focus on this kinetic level model and refer to it as has been done in the current literature [Hill as the AK 
model. It may be regarded from the ecological viewpoint merely as the familiar SI model extended to include spatial 
resolution and diffusive transport. From the mathematical point of view it may be said to represent a system obeying 
the Fisher equation [lol | with internal states representing infection or its absence, respectively. While near-trivial to 
conceive, this model has had considerable success in the short time since it was proposed for the Hantavirus It 
has led to qualitative and semi-quantitative success in explaining observations such as spatio-temporal patterns in the 
epidemics. These patterns are associated with correlations between periods of precipitation and epidemic outbreaks, 
and with the spatial location of refugia-regions of the landscape in which infection persists during off-periods of the 
epidemic '5, ^ . Other applications of the model include the detailed understanding and control of traveling waves of 
infection , fluctuations arising from the finiteness of the numbers and discreteness of the population of the rodents 
[ill HI , environmental effects , curious switchirig effects that have been predicted to occur 14| , and extensions to 
unrelated systems such as bacteria in Petri dishes [isl [l6l | . 

The success of and interest in this so-called AK model represented by Q led to recent attempts to devise practical 
prescriptions for the extraction of the parameters constituting the model from field measurements. The mouse birth 
rate b and the death rate c are obtained from field observations without too much trouble and are generally considered 
to be constants in space and time. With some effort, reasonable estimates of the environment resource parameter 
K {x,t) as a function of location and time can be obtained by counting food (such as nuts and water) available to the 
mice in the different locations, as well as by acquiring aerial photographs of the vegetation cover. Relative, rather than 
absolute, quantification of K is possible in this way. Observational collection of data concerning the aggression rate a, 
through which infection is thought to be transmitted during mouse-mouse encounters, turns out to be so difficult that, 
at least at the present moment, it must be considered an adjustable parameter. It was, however, possible to focus on 
the imp ortant parameter, the mouse diffusion constant D, and to obtain it quantitatively from field measurements 

[il[il[i3. 

The basic idea behind the extraction of D was to regard to the extent possible that the mouse movement is a simple 
random walk, and to extract D from records of the movement through the use of the well-known proportionality 
of the mean square displacement to Dt. The details of the theory [l9| and the iniplementation of the prescriptions 
obtained from the theory to mark-recapture observations carried out in Panama |l7j and New Mexico |3| may be 
found in our recent work. The important and perhaps surprising conclusion that emerged from that work was that 
the mouse mean square displacement, which grows linearly with t for short times, is found to saturate at large times. 
The appearance of a length scale in the random movements of the mice may be ascribed to the fact that animals 
typically move near fixed locations (burrows) for reasons of shelter and security [13, [^ [ll, but it could also be 
ascribed to the fact that mark- recapture observations employ a limited region of space where the traps are laid out. It 
is possible to show analytically 0| that either of these factors could independently lead to the saturation of the mean 
square displacement. A disentangling of the two length scales was possible and led to explicit deductions of both the 
diffusion constant of the mice and their home ranges. Thus, reasonable realistic extracted values of the home range 
size L for different types of mice in different environments were found to be between 50 and 120 m for Zygodontomys 
brevicauda in Panama and about 100 ± 25 m for Peromyscus maniculatus in New Mexico. The respective values of 
the diffusion constant D of the mice turned out to be 200 ± 50 and 470 ± 50 meters squared per day. 



3 



This quantitative information has provided the impetus (indeed, necessity) to generaUze the AK model expressed 
in Eq. ^ to incorporate home ranges. Kenkre has suggested several different generalizations for this purpose. 
One simple way of incorporating home ranges in our model of epidemic spread is to add potential terms to Eq. Q . 
Such an analysis, carried out by Maclnnis et al. ^^23] has resulted in modifications in the AK predictions for refugia 
sizes and shapes. Another, simpler, model modification in the AK equations has led us to apply the so-called 
MontroU defect technique 0| to deduce memory-possessing variations of the AK equations on the one hand and 
time-dependent diffusion constant variations on the other. Perhaps the most fertile model that has emerged from the 
work on the determination of motion and demographic characteristics is the one that we discuss in the present paper. 
It is set out in Sec. ^] below and an examination of its consequences form the rest of the paper. 



II. A GENERALIZED MODEL FOR THE SPREAD OF THE HANTAVIRUS 



Let us consider the dynamics of two types of mice, stationary and itinerant (and susceptible and infected in each 
category). The stationary mice are the adults that move within their home ranges and do not stray far from the 
burrow. We have termed them 'localized adults' in the title of the present paper. For the sake of simplicity as well 
as to emphasize the new features that appear as a result of confinement, we neglect here the fact that home ranges 
may overlap and lead thereby to transmission of infection. The itinerant mice are the subadults (called 'itinerant 
juveniles' in the title) that must leave to find their own home ranges. Adults do not move because their home ranges 
are considered of negligible extent for the purposes of this description. They die at the rate c and have the standard 
logistic competition interactions with the environment controlled by the environment parameter K. They may be 
infected or not, the only possibility of their contracting infection being when an infected juvenile visits their home 
range. If infected, they may transmit infection to a susceptible juvenile if it visits their home range. Adults are not 
born but juveniles turn into adults. This happens when a juvenile finds an appropriate site to settle down in, which 
then becomes its burrow. 

Juveniles are born at a rate b from the adult population. They are mobile, their motion being diffusive. They may 
acquire or transmit infection to other adults or other juveniles on encounter. If they find an appropriate site they 
turn into adults and become immobile as described above. They also have environment competition rates with the 
rest of the mice. In order to focus attention on special features of our generalization, we neglect the death rate cb 
of the juveniles and allow their population to be depleted only through the competition term and their conversion 
to adults through growth. In a companion paper |27| we have included the death rate of the juveniles and indeed 
examined the specific effect of changes in that rate induced by the existence of predators in the open field. 

The characteristics that we have described above suggest that the AK model in be replaced by 

at K[x,t) 

^^f/^ = bA-CBB,-^^^±^-aBM^^B,)+D^^B^-G{x)B,, 
at K\x,t) 

= -cAi — — h aAsBi + G{x)Bi, 

at K(x,t) 

dA,{x,t) A,{A + B) Aj^^rt-M^ (o\ 

= -cAs — — aAsBi-V G[x)Bs. (2) 

at K[x,t) 

where A and B (without suffixes) denote the total densities of the adult and juvenile mice respectively, the suffixes 
i and s represent infected and susceptible states as earlier, and the last terms in each equation describe the settling 
down of the juveniles into their own homes, accompanied by their conversion into (static) adults. The rate of such 
conversion is G{x). This G{x) is non-zero only if x lies in the 'green pastures', that is the spatial regions that the 
juveniles find suitable as their home ranges. Note that there are no spatial derivatives in the equations for the adults 
because the adults do not move, that being an extreme representation of their confinement to their home ranges. 

The model represented in Eq. Q describes the processes at a kinetic level as does the AK model in Eq. (Q. This 
means that the key quantities are mice densities and that the evolution is described via kinetic equations such as 
the Fisher equation ^ ^^^^ detailed description, with no spatial resolution included, is provided by a mean 

field model in which the key quantities are the total mice numbers J A{x,t)dx and J B{x,t)dx where the integrals 
are over the entire landscape. A more detailed description than in kinetic models is provided by configuration master 
equation approaches (see for example jlEES) which include fluctuation effects. In the present paper we will analyze 
both the mean field and the configuration master equation approaches, reserving the middle-level kinetic treatment 
for a future publication. 
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Our paper is laid out as follows. In Sec. IIIII we present an explicit analysis of the mean field model, and the 
simulation treatment of the configuration level approach. In Sec. IIVI we present essential results and their discussion. 
In Sec. El we present concluding remarks. 

III. ANALYSIS: MEAN FIELD AND SIMULATION TREATMENTS 

In this section we treat the model of Eq. ||2J) first by simplifying it to the mean field level of description and then 
by augmenting it to the configuration level description. For simplicity we keep K time and space independent. 

A. Mean field Description 

The mean field description focuses on the time evolution of the integrals of the densities in a kinetic description 
such as that of Eq. (01 and thereby loses space resolution. Called in some contexts the 'well-stirred limit', the mean 
field description may be considered as the limit of Eq. JSJ for an infinitely large diffusion constant D. In passing 
from the kinetic to the mean field description, the single infection quantity a in Eq. (|2Jl results in two corresponding 
quantities: ao and ai. The former (latter) refers to the transfer of infection between an adult (juvenile) and a juvenile. 
Both quantities result from the combination of the infection event and the motion process. The latter is explicit in 
Eq. ||2Jl but not in Eq. ||2Jl. The difference between ag and Oi is precisely the difference between the expression AttRD 
and SttRD which describes the capture and mutual annihilation rates respectively in the literature on excitons in 
molecular crystals 3Q_, Although the precise relationship of oq and ai would depend on the relative importance 
of the mice motion and infection processes, we will take oi — 2ap for simplicity in the rest of the paper in keeping 
with the extreme limits considered in other literature contexts '3.^ . 

The juveniles transmit infection among themselves with rate per unit density ai, get infected and transmit infection 
to the adults with rate per unit density ao, struggle for resources all over space through the environment parameter 
K and become adult with a growth rate g. This growth rate, introduced to represent a juvenile settling into an 
unoccupied home range and growing into an adult, is proportional to G{x) of Eq. Q. The other terms can be 
interpreted as explained above for Eq. ^ . The coupled set of equations for normalized quantities in the mean field 
description is 



dr 

dr 
dAi 

dAs 



-7S, + aoBsA, + aiB.Bs - B., [A + B), 

^jBs + (3{As+ A^) - aoBsA^ - aiB.Bs -Bs{A + B), 

-Ai + aoAsBi + -fB, - A, (A + B), 

-As - aoAsB^ + jBs -As{A + B). (3) 



where each script character denotes the ratio of the quantity described by the corresponding Roman character and cK . 
Thus, Ai — Ai/{cK). We also write t = ct and the dimensionless parameters are now 7 — g/c, (3 — b/c^ ao — Kao, 
and ai — Kai. As explained above we take ai — 2ao. 

From Eq. ||2Jl it is easy to recognize that the total adult and juvenile populations, respectively, A = Ai + As and 
B = Bi -\- Bs, obey the following evolution 

= (3A--fB-BiA + B), 

dr 
dA 

= -A + -iB-A{A + B). (4) 

dr 

Even if A and B are both limited by a quadratic saturation, it is evident from (0J that the sum A-\- B does not obey 
the standard logistic equation. This fact, as we will see in Sec. IIVI gives some qualitative differences of the steady 
state parameters dependence of (O from the AK model. The steady state values A and B are given for /? > 1 by 
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with ^ = 2(/3 — 1)7/(1 + 7)^. The situation /3 > 1 represents the juvenile birth rate b being larger than the adult 
death rate c. For the opposite situation, /3 < 1, when the adults die quicker than the juveniles are born, the trivial 
solution A = B ~ emerges from Eq. (@}. These two steady states exchange their stability as /3 crosses the value 1 
clearly indicating the presence of a transcritical bifurcation at /? = 1 from zero to non zero population density. 





FIG. 1: Effect of variation in the growth rate on mice densities as given by the mean field description. Steady state mice 
densities are plotted in (a) against the normalized growth rate 7 (= g/c) at which juveniles grow into adults, and the evolution 
of the total density of infected mice, normalized to its initial value, is plotted in (b) against the normalized time ct. Both 
exhibit non-monotonic behavior, for instance, a rise in infection as 7 is increased from zero and a decay beyond a certain value. 
System parameters have been chosen to be ao = 1.1 and P — 15. In (a), the four mice densities have been shown in the main 
figure and the total infected and susceptible densities, sums of adult and juvenile contributions, are shown in the inset with 
a:xes identical to the main figure. In (b), the ten curves are for 7 incremented by 1 from to 9 and go down the graph in 
order for short times but for long times approach equilibrium values that are not in the same order. Notice also the peculiar 
non-monotonicity in time. 



With the help of the non-zero steady state values A and B in Eq. Q it is possible to obtain (see Appendix) the 
steady state solutions of Eq. |j3Jl. The system has the trivial solution {Ai = As = Bi = Bs = 0) when /3 < 1, and 
two possible steady states for /3 > 1. They represent respectively the non-infected and infected phase. The former is 
given by 

B^ = 0, 
Bs = B, 
A, = 0, 

As = A, (6) 



while the latter is given by 



B., = 



+ 8£^T 



4ao 

Bs = B-B,, 

B, (7 + aoA) 



a^B. + l+A + B' 

wherein = 2(7 + 1) + ^4 (3 + ao) + B (3 - 2ao) and f = a^B (7 + a^A) - [7 + ^4 + F (1 - 2ao)] {l+A + B). A 
stability analysis (see Appendix) shows that the solution sets represented by Eqs. Q and (jT)) exchange their stability 
through a transcritical bifurcation when the normalized infection rate is larger than a critical value Uc- By determining 
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when £" = we can calculate the critical value as shown in the Appendix. Equivalently, by studying the non-normalized 
mean field equations, we can calculate the critical environment parameter Kc'- 



j + 2{l + A + B) 
2aoA 



^ AA {-i + A + B) {l + A + B) _ ^ 
B[7 + 2(l+:4 + 6)]' 



(8) 



Here A~\- B is the total population at steady state, i.e., the carrying capacity of the system normalized to cK: 



A + B = {P-l) 



7 1 

1 + 7 i 



(9) 



An infected phase exists for K > 



B. Simulation Description 



It is well-known |33j that spatial aspects, not accessible to mean field theory, are (obviously) very important in 
epidemiology, generally in ecology, and that they require a description, capable of addressing spatial correlations. One 
way to treat spatial resolution is to adopt kinetic level approaches as in the AK model, while another, a more detailed 
way, is to adopt an approach based on the evolution of the full configuration states. We have elected to choose the 
latter in the present paper. Analytic solutions are typically impossible in such an approach except for oversimplified 
models. Hence we resort to simulations as done previously [Ol for the AK model. 




FIG. 2: Effect of variation in the growth probability as given by the simulation description. Compare with Fig. Q Steady 
state mice densities are plotted in (a) against Pg and the time evolution is shown in (b) where the ten curves correspond to 
increments of 0.1 in Pg. Contrary to the mean-field model (see Fig. 0, where infected population as function of the growth rate 
7 always displays a non- monotonia behavior, here the infected population has a monotonia increase with Pg. The parameters 
for both (a) and (b) are Pc = 0.01, Pb = 0.011, Pa = 0.3 and Pk = 0.99. 

Our simulations are carried out on a L x L square lattice with each site of the lattice corresponding to a small 
region in the landscape. Moderately large lattices (with a total of 2^** sites) have been used in the simulations. The 
four subclasses of mice, (adults and juveniles in susceptible and infected states) change their numbers, as time evolves, 
in accordance with rules which represent the model under consideration. At each time step, the juveniles may move 
but the adults not, the probability for the diffusive (random walk) motion being 0.125 for any of the eight directions 
of the square lattice. We consider the time step scaled to the diffusion constant in this manner. An adult, infected 
or susceptible, gives birth to a susceptible juvenile with probability P^. An adult dies by aging with probability Pe- 
lf two or more mice meet at a site, one of them may die with probability 1 — Pk- If a susceptible mouse occupies 
the same site as an infected mouse, the former has probability Pa of getting infected in the next time step. And if a 
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juvenile mouse finds itself at a site without an adult, it grows up and settles at that site with probability Pg. These 
rules represent a simplified version of Eq. ^ augmented to the configuration level. It is simplified in that the regions 
('green pastures') where the juveniles may settle down have not been marked but the process has been represented 
through a probability of conversion. We have carried out full scale simulations which take into account the spatial 
extent of the green pasture regions and will discuss the results elsewhere. 

The simulation description thus has the parameters, Pt, Pc, Pg, Pa, and Pk, in correspondence to the respective 
rates b, c, g, ag, and K of the mean field equations analyzed above. However, the correspondence is not straightforward 
in all cases, as expected. Extensive computer simulations were performed. We found that the system reaches a steady 
state, after a transient. The main quantities we analyzed were the densities of the four populations of mice and the 
existence or not of infection in the steady-state. Densities are defined as the total number of mice in the lattice divided 
by the total volume (or area) of the lattice. Since the existence of more than one adult per site is not allowed, the 
density of adults lies between and 1. By contrast, the juvenile density may exceed 1. 



IV. RESULTS 



The main difference between the mean field and the simulation results lies in the effects of the quantities that govern 
the growth of juveniles into adults, rate g and probability Pg, respectively. We show this dependence respectively in 
Figs, n and 13 From Fig. 1(a) we see a non-monotonic dependence of the steady state infected population densities 



1(b) where the time evolution for the entire infected 



as 7 increases. This non-monotonicity is also evident in Fig. 

population is shown for different values of 7. In the limit 7 — > in Eq. ^ A, B —^ since no new adult can be 
'born' (i.e., produced by conversion of a juvenile through growth), the juvenile population cannnot be regenerated 
and dies out because of competition, and eventually no population can be sustained. The juveniles are responsible for 
spreading the infection to the adults and a larger number of them tends to increase the infected population. However 
if the juveniles convert into adults too fast, there are less mobile carriers of infection, the number of infected juveniles 
decreases (note that juveniles are born susceptible and never infected, this being a Hantavirus characteristic), and 
eventually the infection disappears. In other words we can say that as 7 grows, the critical environment parameter 
Kc eventually becomes smaller than the system K. 
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FIG. 3: Steady state densities versus the environment parameter plotted from the mean field description in (a) and the 
simulation description in (b). The x-axis coordinate in the inset and the main figure is K in (a) and Pk in (b). Other 
parameters are taken to he P — 1.5, 7 = 1 and ao = 1 in arbitrary units in (a) and Pc = 0.01, Pb ~ 0.014, Pg — 0.1 and 
Pa ~ 0.3 in (b). The behavior is similar in both levels of description. When the environment parameter is large enough, an 
infected phase emerges. Once that happens, the increase of the population is due only to a larger density of infected animals, 
the susceptible population decreasing its overall density to a constant. The susceptible population in the simulation description 
becomes zero for Pk smaller than a critical value that depends primarily on the value of Pb relative to Pc- If Pb is sufficiently 
close to Pc, as is the case shown here, competition and adult death processes drive the entire system to extinction at a nonzero 
value of Pk- 
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In Fig. |2(a)| we show the Pg d epend ence of the steady state populations in the simulation description. This is 
a noteworthy difference from Fig. |l(a)| a monotonic increase of the infected population. This difference has to be 
ascribed to a spatial correlation effect. Since there can be only one adult per site, an increase in Pg has also the effect 
of reducing the number of available sites for the juveniles. For sufficiently large values of Pg, clusters of adults start to 
form, creating confined regions in the landscape where the juveniles are constrained to roam. This in turn increases 
dramatically the average time necessary for a juvenile to find an available site to settle down. Unfavorable effects of 
the growth of juveniles into adults observed in the mean field theory do not therefore occur in the simulations as a 
result of the spatial correlations set up by the cluster formation in the adults. The monotonic increase of the infected 



population as function of Pg can also be observed in Fig. 2(b) where the time evolution for different Pg values is 
depicted. In order to verify the validity of the above explanation of the difference in the mean field and simulation 
predictions, we studied a pseudo-model, intermediate between the two descriptions. In the pseudo-model, a juvenile 
can move to an arbitrary position and not only to a nearest-neighbor site, thus being able to jump the barriers set 
up by the clusters formed in the adult population. Careful simulations we have carried out show that, indeed, the 
steady state populations have an infected phase that decays to zero beyond a critical value of Pgi the pseudo- model 
predicts the same qualitative behavior as the mean field theory. 

The struggle for resources described at the mean field level by the environment parameter K is here represented by 
the probability of survival Pk- the probability of dying via competition for resources is 1 — Pk- In Fig. 3(a)|we show 



the steady state populations as function of K in the mean field description and in Fig. 3(b) the corresponding results 
as function of Pk) in the simulation description. The behaviour in the two cases is similar: beyond the critical value 
of K (Pk) the infected population increases with a sublinear dependence on K (Pk)- 

V. CONCLUDING REMARKS 

From a perspective of concepts that have been successful in physics, our present model of the mice assembly may 
be described by the term liquid-solid because it describes a class of mice that move freely as do the molecules of a 
liquid, and another class of mice that move in the neighborhood of fixed positions (the burrows) as do the molecules 
of a solid, which vibrate around fixed lattice sites. Therefore, for the purposes of the following discussion, we will call 
our model the liquid-solid (LS) model. 

Although the LS model was constructed by generalizing the structure of the AK model, the two have really only 
one parameter truly in common: the environment parameter K. The birth rate b in the AK model describes the 
emergence of new mice from all mice whereas it produces only juveniles from adults in the LS model. The death rate 
c in the AK model is similarly applicable to all mice but in the present version of the LS model it applies only to the 
adults. In addition to depletion via competition for resources, the disappearance of the juveniles is assumed to occur 
only through the growth rate g when they grow up into adults. The diffusion constant D describes the motion of all 
mice in the AK model but only of the juveniles in the LS model, the adults being stationary. Therefore, we make 
comparison comments regarding the two models by looking at how the infection depends on K. For simplicity we 
discuss this comparison at the mean field level. 

At first sight, it might appear that the LS model does not reduce to the AK counterpart in any situation. However 
there exists one such limit. Let K be infinite so that the carrying capacity is infinite also. To eliminate runaway 
solutions in the steady state, let us take the birth and death rates equal in the AK model so that Bak — cak, and 
let us correspondingly take the birth-death ratio /3 = 1 in the LS model. Consider now the limit in which the growth 
rate g of the juveniles and the death rate c of the adults in the LS model are much larger than the rate of infection 
between adults and juveniles but such that ^ — g/c is finite. In such a situation the adults constitute an isolated 
reservoir of the infection: they do not play any role in spreading it since an adult is infected only if it 'grew up' from 

1/2 

an infected juvenile. In particular, it can be shown that by considering ai = qak and 7 = \{1 + icAK) ' + l \ /2, 
where qak is the infection parameter in the AK model, our present model at steady state gives exactly the same 
analytic dependence of the AK model for the infected and susceptible population as function of cak and uak- 

This equivalence suggests a simple way to compare the two models, when K is finite, by considering that no 
transmission of infection occurs between adults and juveniles, i.e., oq is taken to be zero in the LS model. In such 
a scenario, the two models have the same qualitative dependence of the infected and susceptible populations as 
function of the environment parameter: for K > Kc the susceptible population remains constant while the infected 
population increases linearly. In other words, by increasing K, the additional population (proportional to K) that 
the environment can sustain eventually becomes infected at steady state. In Fig. 0]we show this dependence for the 
AK model at the mean field (Fig. |4(a)|) and at the simulation description (Fig. 4(b) I, respectively. The behavior is 



similar in both levels of description. When the environment parameter is large enough, an infected phase emerges. 
Once that happens, the increase of the population is due only to a larger density of infected animals, the susceptible 
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population remaining constant. By contrast, as previously shown in Fig. the infected populations increase while 
the suscpetible populations decrease to a non zero value as K > K^. At the mean field description it is also evident 
that the infected population increases only sub-linearly as function of K and not linearly as in the AK model. 

We indicate how our theoretical predictions may be used in conjunction with observations in the field to test the 
validity/ applicability of the LS model. The qualitative differences between the AK and LS models discussed above 
could be exploited in analyzing data from mark-recapture observations. In such observations, traps are set up in regions 
in the landscape and mice are caught, examined, marked, and released. Information about a variety of features is 
gathered including infection status and age. Juveniles have clear physionomical characteristics that distinguish them 
from the adults |2j| . If all other effects are controlled, supplementing additional food homogeneously over the terrain 
would be a way to increase the environment parameter K. Experiments that exploit this feature may allow one to 
determine if, in the presence of infection, the susceptible population remains constant or decreases as the amount of 
food is increased. This would allow us to establish whether the juveniles are indeed the main carrier of the disease 
and whether augmenting the AK considerations to the LS model is the correct way of analyzing the spread of the 
Hantavirus infection. 
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FIG. 4: Steady state densities of the Ak model as function of the environment parameter plotted from the mean field description 
in (a) and the simulation description in (b). The x-axis is the environment parameter K in (a) and the probability Pk in (b). 
Other parameters are taken to be /3 = 1.1, and aAK = 30 in arbitrary units in (a) and Pc = 0.01, Pb = 0.02, and Pa = 0.6 in 
(b). 



Since in both models the value of Kc depends on all the other parameters, there are obviously situations in which 
one model predicts infection while the other does not. However, if we compare the amount of infection as K increases 
beyond our analysis shows that the AK model gives more infection than the LS model. Surely, this is to be 
expected since, in the latter, part of the population (adults) is stationary and transmit infection less efficiently. This 
is an important consequence of the existence of home ranges determined quantitatively in our recent work p!7lll8l[l9l |. 
Our analysis allows us to quantify that consequence, i.e. to determine how much the reduced motion of the adults 
decreases the transmission of infection compared to the AK model. While, for simplicity, we have considered here the 
extreme limit of zero overlap between neighboring home ranges, to what extent the degree of home range overlap will 
change our conclusions is an open problem and the subject of our current investigations. 

To keep our analysis focussed on essentials, we have considered in the present paper the case in which the juvenile 
population is depleted only through growth into adults and not through death. We have carried out mean field as 
well as spatially resolved studies of the situation when this simplification does not apply. This is important because, 
in their itinerant attempts to find suitable places for their own home ranges, juveniles are surely exposed to predators 
that kill them. Details of the investigation will be given elsewhere but we report here that we have uncovered 
the epidemiologically noteworthy possibility of buffering the transmission of the infection by introducing predators in 
the landscape. Sustained by a sufficiently large K (environment resources), a large susceptible population of mice can 
exist but with no infection as a consequence of interaction with predators. 

Among additional avenues of theoretical research based on the LS model are the study of spatial correlations based 
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on modern techniques [SJI, the extension of our analysis by specifying the spatial dependence of rates G{x) in Eq. (0) 
so that the effects of the location of the 'green pastures' can be ascertained, and the comparison of our predictions 
to observations in the field concerning infection spread. The latter effort is particularly important because it will 
allow us to explore the applicability and practical relevance of home range inclusion in the theory of the spread of the 
Hantavirus. Research in all these directions is under way. 
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VI. APPENDIX 



Since the total adult and juvenile population in Q is known in terms of the parameters ^ and 7, the set of equations 
(O at steady state can be simplified considerably to 

aoB,:A^ + alBs^^-B^{"/ + A + B} = 0, 

B, + Bs = B, 

(7 + aoX)B, - A(l+^ + 6) = 0, 

A + As = A. (10) 

The system Hl()(l has only three possible solutions with the four variables larger or equal to zero. The trivial one 
implies A = B = 0. The other two solutions can be obtained by reducing through substitution Eq. H10|) to the 
following polynomial equation in Bi 

Q^o + "0 + ai) + B[ao + ai (1 — ao)] + 2^ao + ai} 
+ [7+:4 + B(l-ai)] {l + A + B) - aoB {-/ + aoA)} -0, (11) 

from which it is easy to obtain the solutions shown in Eq. © and {Tj) . The study of the sign of Eq. 111|) gives the 
condition for the existence of an infected phase 



+ 2(1 + A + B) 

Oic = — = 

2A 



^ ^ 4A{-f + A + B) (l + A + B) ^ 
\i B[-f + 2(1 + A + B)]^ 



(12) 



which can be converted to a Kc as written in Eq. JSJ. Notice that the third root of the polynomial in 1)11(1 can be 
shown through a numerical study to be always negative and it is thus discarded. 

The stability analysis of the three solutions of ^ is done by calculating the Jacobian of the system at steady state 
J (Ai, Asj Bi, Bg) 



J = 



/ -(l + A + B + A,) _aoBi~A 

-A _ - {l+A + B + As+ aoB^) 
aaBs - Bi_ -Bj_ 
V f3-{ao + l)Bs f3~Bs 

aoAs - Ai_ -A_ \ 

_ -(ao+J_)A _ _ ^-^'^ _ 

2aoBs- {j + A+_B + B^) jvqA, + (2ao - l)Si 

-{2ao + l)Bs -[-f + A + B + Bs+ao{A^ + 2B^)] J 



(13) 
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The trivial solution has the following four eigenvalues 



Ai = -1, 

A2 -7, 

7 + 1 



A3 = 
A4 = 



2 

7+1 



+ 2e + 1 

(v/l + 2e-l 



(14) 



from which it is evident that A4 < if /? < 1, while all the other eigenvalues are always negative. The trivial solution 
is thus stable if /3 < 1 and it becomes unstable when f3 > 1. 

The polynomial characteristic P{X) of the Jacobian H13|l associated with the solution is equal to the product 
Pi(A)P2(A) where 



Pi(A) = A2 + A[1 + 3(^ + 6) +7] 



(1 + 2^ + B) {A + 2B + j) - (B- P) {A - 7) 



and 



P2(A) = + X[l + 2{A + B) -2aQB + j] 

+ {1 + A + B) {A + B-2aoB + -f) ~ alAB. 

The eigenvalues associated to Pi (A) are given by 



AlH 



1+7 



1-3^1 + 2^+^2 (l+e+ VTT2I 



(15) 



(16) 



(17) 



with Ai_ < for any /3 and 7 and Ai_^ > when /3 < 1. The eigenvalues associated to P2(A) are negative for ao < ol^ 
and they become positive when ao > etc- The solution defined in Eq. © is thus unstable if either /3 < 1 or if ao > Ofc- 
The study of the sign of the eigenvalues associated with the third possible steady state is done numerically and it 
is possible to show that Eq. represents a stable steady state if /3 > 1 and ao > etc and become unstable if either 
/3 < 1 or ao < ac. The mean field description from the dynamical point of view has thus two transcritical bifurcations: 
one when the growth rate h equals the death rate c and the other one when the infection rate ao equals the critical 
infection rate or similarly when the environment parameter K equals the critical environment parameter K,.. 
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